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ABSTRACT 

A  general  analysis  of  coplanar  lines  on  a  dielectric  substrate  is 
presented.     The  field  equations  are  written  in  terms  of  Hertzian  poten. 
tials  and  boundary  conditions  are  applied  to  generate  a  system  of 
linear  equations.     These  equations  are  then  Fourier  transformed  and 
reduced  to  a  form  where  a  computationally  efficient  solution  can  be 
obtained  by  the  method  of  moments.     Any  desired  accuracy  may  be 
obtained  by  expanding  the  transform  of  the  current  if  the  problem  is 
one  involving  strip  condxictors  or  the  electric  field  if  there  are  slots. 
A  one  term  expansion  produces  very  good  results.     Both  wavelength 
and  characteristic  irapedance  of  the  transmission  line  structure  are 
obtained.      Theoretical  and  experimental  results  are  presented  for 
coplanar  strips. 
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I.     INTRODUCTION 

During  the  past  several  years  a  variety  of  coplanar  microwave 
transmission  lines  has  been  developed  covering  a  wide  spectrum  of 
applications.  However,  there  still  remained  several  possible  struc- 
tures to  analyze,  such  as,  coupled- slot  lines,  coplanar  strips,  etc. 
Also,  theoretical  analytic  methods  more  accurate  were  sought  to 
reduce  the  errors  obtained  by  the  more  classical  methods,  such  as 
the  electrostatic  approximation,    the  niethod  of  partial  iinages,    etc. 

The  present  work  is  concerned  with  two  parallel,    coplanar  strips 
over  a  dielectric  substrate,    the  other  side  being  bare,    as  illustrated 
in  Fig.    1.     The  structure  is  uniform  and  infinite  in  both  x  and  z  direc- 
tions; it  is  also  assumed  that  the  substrate  material  is  lossless  and 
its  relative  permittivity  is   C 

For  coplanar  strips  (GPS)  to  be  practical  as  a  transmission  line, 
radiation  must  be  minimized.     This  is  accomplished  through  the  use 
of  a  high  permittivity  substrate,    which  causes  the  coplanar  strip 
wavelength  to  be  small  compared  to  the  free- space  wavelength,    and 
thereby  results  in  the  fields  being  closely  confined  to  the  strips  with 
negligible  radiation  loss. 

A  spectral  domain  transform  method  was  suggested  by  Professors 
Itoh  and  Mittra  [Refs.    1,    2]  which  yields  an  exact  solution;  this  method 
together  with  the  method  of  moments  [Ref.    3],    is  used  in  this  present 
work. 
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Fig.    1.     Parallel  coplanar  strips  configuration 


II.     DISPERSION  CHARACTERISTICS 

A.      FIELD  AND  BOUNDARY  CONDITIONS 

Let  the  following  electric  and  magnetic  fields  exist  along  the 
propagation  or  z_direction,    as  justified  by  the  Hertzian  vector  poten- 
tial functions  described  in  Appendix  A, 

^        ^=    ^  ^  (Al) 

4i^    :=    L'    $^'U^  ;  (A2) 

where    ^     is  the  propagation  constant,    and  since  it  was  stated  that  the 
substrate  material  was  lossless,    it  follows  that 


^=J(^     . 


(A3) 


From  the  field  expressions  in  equations  (Al)  and  (A2),    all  other 
components  of  electric  and  magnetic  fields  can  be  derived  from 
Maxwell's  curl  equations,    as  shown  in  Appendix  B,    as, 

(A4) 

(A5) 

(A6) 

(A7) 


^y      )    2 

t"  1  ^' 

J/       (  y  ^"''      • 

r  \ ." 

Applying  boundary  conditions  at  the  interface  between  regions  2 
and  3,    tangential  electric  and  magnetic  fields  must  be  continuous; 
therefore,   at  y  =  0, 


(A8) 
(A9) 
(AlO) 
(All) 


Similarly,    at  the  interface  between  regions   1  and  2,    tangential 
electric  fields  must  be  continuous  and  tangential  magnetic  fields 
discontinuous  by  corresponding  surface  current  densities;  therefore, 
at  y  =  d, 


fei,  (x,a,a)  =    6^J,(t,c*,^) 


[   o  dseuhcre 


y^ 


^x^(x,a,-ij~  Ux^  (^,ci,e)  = 


3^^A  ^  on  strljps 


(A12) 
(A13) 

(A14) 

(A15) 


Also,    the  electric  fields  will  exist  only  in  the  dielectric  part  of 
the  interface  and  can  be  expressed  as, 
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^,(r,,a,i)^{  ^  (A16) 

.  <^  on  strifes 

Substituting  the  field  expressions  of  equations  (Al)  through  (A7) 
into  the  boundary  conditions  expressions  of  equations  (A8)  through 
(A17)  one  obtains, 

(A18) 
.xC^l  .  itel  5-(k)  _(ffl 


V  ^  C^c.)  -j«);).  ^±  (..d)  .  r  ^  (.,d)  -  ju>^.   f^  (x,d^  (A23) 


W  ^?^  (-^.cJ)  -  tcj^  $^i^'  (%,d)    =  J^  (^^  (A24) 

7  (A25) 

ic?^'*^  (^,d^   =  e^M  (A26) 


(A27) 
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B.      SPECTRAL  DOMAIN  TRANSFORM 

All  potential  functions  must  satisfy  the  following  relation: 

(Bl) 
where, 

'  (B2) 

One  introduces  the  Fourier  transform  to  the     oc  -domain,    as 
suggested  by  Itoh  and  Mittra  [Refs.    1,    2],    via, 

Therefore,    one  can  find  the      csc -domain  Fourier  transform  of 
equation  (Bl)  as, 

(B4) 


(B5) 


Furthermore,    one  can  define  propagation  constants  for  each 
region. 


7. 


(B6) 
where, 
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From  planar  dielectric  substrate  theory,    the  effective  permittivity 
is  related  to    €       and     €  j  as, 
€.  <6^^<€2  (B7a) 

Therefore,    analyzing  the  expressions  for  propagation  constants 
in  equation  (B6),    for  region  1  or  3, 

where     A     is  the  dielectric's  effective  wavelength  and  is  related  to 
the  free- space  wavelength     A      as, 

^  <  >^'  <  >^   .  (B8a) 

Replacing  this  last  expression  into  equation  (B8),    one  can  find 
lower  and  upper  bound  values  for    0    i   ^^    ^    ^  as, 

o<^+(€,.l)(-^\^  ^y>^^2  (B8b) 


therefore,     0        ^-^d    0       are  always  real  quantities,    independent  of 
the  values  of      0<    . 

Similarly,    for  region  2, 


N-«H,^-t^.^.(f;-(^y 


(B9) 
where,    as  before,    the  upper  and  lower  bound  values  for    Cx:    are, 

\    ^  '    .  (B9a) 

It  is  clear  that   Y  _  will  be  imaginary  for  sinall  values  of    oc    and 
real  for  large  values  of     o«i"    .     Specifically,    %       will  be  imaginary 
for 
«<  <    ■—  V^r-i  (B9b) 
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and  real  for 

0<>  ^  \j^-l      ^  (B9c) 

Solving  the  potential  functions  differential  equation  (B5)  for  each 

of  the  three  regions  one  obtains,    as  derived  in  Appendix  C,  -  Electric 
field  potential  functions: 

X,<«)              ,          ^^t)           -^-r^r^^  (BIO) 

<I>^*'(»c,.^)  «  B/*'(6<)  &/n  r^^  +   C^*^'(»c)  cos  YtU  ^   Yj-cStu.  (Blla) 

$3  (\'^)=  B   (K)  e   ?  (B12) 

-  Magnetic  field  potential  functions: 

* (fc)  ,         .        .(U^          -riff-d)  (B13) 

$  J.  (^ty)  -    B,  (et)  s/o    »V^  +   C^*"'^**)  COS  y^if  ,  Vj  =  <Qut  (B14a) 


where,    clea: 

rly 

v;- 

Va 

L, 

=  tc. 

and. 

(B16) 
(B16a; 
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^'•^.^^;p!'*^^.^  (B17) 


(B17a) 


where  both,     B     g^y,    and  B^  ,    are  evaluated  at    V        =  j  Y  ^ 

In  all  the  following  discussion,    equations  subscripted  with  lower 
case  "a"  will  refer  to    0        being  imaginary. 

The  set  of  equations   (BIO)  through  (B15)  can  be  interpreted 
physically  as  follows: 

Electric  and  magnetic  fields  are  more  concentrated  in  the 
dielectric  substrate,    region  2,    and  decay  exponentially  toward 
zero  outwards  into  the  surrounding  media,    regions    1  and  3. 
In  the  dielectric  substrate  itself,    electric  and  magnetic  fields 
and,   hence,    energy,    are  more  concentrated  in  the  lower  values 
of      O^    .     Although  there  is  no  direct  comparison  of  the  domains 

Cx^    and    X    ,    a  set  of  plots  of  the  potential  functions  indicate 
that  energy  is  also  confined  to  the  vicinity  of  the  strips. 
An  approximate  graphical  representation  of  the  energy  distribu- 
tion within  the  dielectric  substrate  for  a  fixed  value  of  dimension  y 
is  shown  in  Fig.    2. 

It  is  convenient  to  obtain  the  Fourier  transform  of  equations 
(A18)  through  (A26)  as  follows. 
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Fig.    2.      Energy  distribution  within  the  dielectric  substrate 
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(B18) 


-j..v-§f\»<,oUj«>^2.-^'('',°)  =-J«<ir#»"<''.°)+j«"^-%-<'''°>  (B21) 

|^,'^'^(«,a)  =  t|4'"K=»)  ^         (B22) 

J       -'        '        "J    /         ai{  "  °i (B23) 


«l  1  (B25) 


4;i>/"('''^^  =  ^«^^^  j  (B26) 

'1 


.i.Ya/^^.,a^  -jco^.-^(K,dl  «=  ^J^<)  ^^^^^ 


where  the  derivative  transform  pair, 

^    tl      "  "vi''^  (B27a) 

has  been  applied. 

Upon  substitution  of  the  field  expressions  of  equations   (BIO) 
through  (BIS)  into  the  above  relations,    these  become, 
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(B28) 
ce.)  (U)  te<  (L\  (B29) 

tc?  ^^*^(M  =  ic^^[ a^*'(M  s/oti  rj d  +  c^'^h^,)  cos^  r^al  (bsz) 

(B33) 

+  t*>po^z  \h^'"\tK)  cosh  r^d  +  c'''\k)  s/oh  Yzd] 

(B35) 

I 

icf  A^''(o()  =   i^(6t)  |(B36) 

I 

Equations  (B28)  through  (B31)  can  be  expressed  in  terms  of 

C  and  C       '     ,,    obtaining: 

(oc)  t.-^) 

D'^Vm.(^;)'c'"<.<^  (B39) 
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where,    in  equations  (B40a)  and  (B41a), 


>2=rY, 


Fui'thermore,    one  can  relate  the  constants  C  and  C 

{^)  (c*r ) 

in  terms  of  A     (o^)  and  A       (Oc)  via  equations  (B32)  and  B33), 
obtaining, 


A  relation  between  the  constants  A       {oc  )  and  A       (o^  )  in  terms 
of  the  surface  current  densities   J^C^  )  and  J^{fiC)  can  be  established 
via  equations  (B34)  and  {B35),    obtaining. 


Finally,    applying  equations  (B36)  and  (B37),    the  surface  current 

densities  can  be  related  to  the  electric  fields     ^    {^)  and     ^     (od), 

X  z 

obtaining. 
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(B46) 

where  all  the  intermediate  steps  are  fully  derived  in  Appendix  D. 

C.      PHYSICAL  PARAMETERS 

Up  to  this  point,    the  problem  has  been  stated  in  a  quite  general 
form,    i.  e.  ,    there  is  no  dependence  on  the  actual  physical  configura- 
tion other  than  the  boundary  surfaces. 

In  general,    the  surface  current  densities,    J   (tj^  )  and  J   (Ck:  ),    can 
be  expressed  as  a  train  of  known  basis  functions,    [Ref.    3],    as, 

Jk(«^1-  2  a^  jx-  C«^  (Cl) 

GiCtxW  £   ^^  i^i^^^  (C2) 

Substituting  the  above  expressions  into  equations  (B46)  and  {B47), 
one  obtains, 

Since  the  two  conductors  are  sufficiently  narrow,    one  may  assume 
that  the  surface  current  density  in  the  x-direction  is  zero,    which  will 
modify  equations   (C3)  and  (C4)  as. 
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Furthermore,    although  theoretically  one  should  perform  a  multi- 
term  approximation  to  the  surface  current  density  in  the  z-direction, 
in  this  case  a  one  term  approximation  is  sufficient  since  a  good 
assuinption  of  this  surface  current  density  is  a  square  pulse,    i.  e.  , 

Ml  (c^,^)  ^»  i^, (^"i  =  <ii  ('^^  (C7) 

'  '^  •  (C8) 

or,    substituting  equations  (CI)  and  (C2)  for  the  case  where  i  =  1, 

M,  (.f^,^)  Ji  (tc)   =  g^  (tsci  (C9) 

'  (CIO) 

It  is  clear  that  one  only  needs  to  work  with  either  one  of  the 

above  equations;  consider  equation  (C9),    for  the  present  case. 

As  mentioned  before,    a  good  approximation  to  the  z-direction 

surface  current  density  is  a  square  pulse;  therefore,    it  is  necessary 

to  obtain  the     O^  -domain  Fourier  transform  of  J   (o<  ),    shown  in 

z 

Fig.    3,    as. 


-I-        i 


ctiC 
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J**  L 


jKx 


+  e 

s 

4-^ 

|+*o 

J-f^-J*'*  ->(|-+'^i  J^f  Xi-^vol 


3i[ 


<2  ~  e 


^  6a<{  N  >  -  ^  <io<} 


-<•  <E      -e 


^.(IW;] 


,  Oi   s/o  -I  (s^vJ)  e;o  I  vO 


Therefore,    equation  (C9)  becomes, 


(Cll) 


j<4 


M,  (^,^)  ~   S'V)   ^  (£^-w^  sin    ^  vj    =    <§^(o<:> 


(C12) 


In  general,    for  a  strip  configuration,    the  electric  fields  are  not 

easily  defined;  therefore,    it  is  convenient  to  make  equation  (C12) 

independent  of  vS    (ex).      This  could  be  done  by  taking  an  inner  product 
z 

with  a  function  orthogonal  to    ^    (oc).     Applying  this  concept,    equation 

z 

(C12)  becomes, 


Define  a  suitable  inner  product,    such  as. 


(C13) 


(C14) 
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Fig.    3.      Z  -  directed  surface  current  density 
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A  suitable  weighting  function  W{Cxi  )  could  be  the  complex  conjugate 

of  J   («v  ),    i.  e.  ,    J   (-  o<:). 
z  z 

Also,   by  Parseval's  theorem,    the  right-hand  side  of  equation 
(C13)  becomes  zero  because  of  the  orthogonality  of  the  integrand. 
Therefore,    equation  (C13)  becomes. 

It  is  clear  that  once  equation  {C14)  is  integrated,    the  dependence 
on  the  variable   c<  disappears,    so  one  can  state  that  the  whole  process 
is  really  a  function  of  frequency  and  the  structure's  physical  charac- 
teristics,   namely,    width  of  the  strips  W,    separation  between  the 
strips  S,    thickness  of  the  substrate  D,    and  the  dielectric's  relative 
permittivity    £     • 

Once  the    oc  -domain  Fourier  transform  of  the  surface  current 
density  is  obtained,    expressions  for  the  different  constants,   A  through 
D,    have  to  be  redefined  to  account  for  the  imaginary  character  intro- 
duced by  J   {oc).     It  turns  out  that  all  the  magnetic  potential  constants, 
z 

.(h).       .      ^(h)  (h)  (h)         ...  , 

1.  e.  ,    A       {OC),    B       (O^),    C       (<x),    and  D       (o^),    are  imaginary,    and 

all  the  electric  field  constants,    real.     A  detailed  discussion  of  this 

procedure  is  presented  in  Appendix  E. 

Using  the  expressions  derived  in  Appendix  E,    the  equations  for 

the  various  constants  are  modified  as, 
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a'^Vm-J-— J^l -|J,m|  ,C17) 


c'"(«Wj[|P,1a'''m-PJa«''(m|]  (C19) 


D.     NUMERICAL  INTEGRATION 

Essentially,    the  inain  problem,    as  far  as  numerical  integration 
is  concerned,    is  the  evaluation  of  equation  (C15).     Although  theo- 
retically the  limits  of  integration  are  -  oo    and  +  oo,    truncation  at 
points  corresponding  to  sufficiently  small  values  of  the  integrand  is 
allowed. 

The  numerical  integration  method  used  is  a  modified  Simpson's 
rule.     Consider  a  general  function,    as  shown  in  Fig.    4,    where  one 
divides  the  axis  of  integration    c>i    in  N  segments  of  width    A    ;  the 
ordinate,    corresponding  to  the  half-v/idth  point  of  each  segment,    is 
taken  as  the  average  ordinate  of  each  integration  segment  and  a 
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Fig.    4.     Simpson's  rule  of  nuinerical  integration 
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-i5o         rso 


Fig.    5.     Typical  integrand  plot 
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rectangle  is  formed  with  the  segment's  width.     It  follows  that  the 
total  area  under  the  given  curve  will  be  approximately  equal  to  the 
sum  of  the  areas  comprised  by  the  N  rectangles. 

It  is  clear  that  for  inore  accurate  results  the  width  of  each  seg- 
ment  A    should  be  a  function  of  the  rate  of  variation  of  the  function 
F(o^),    i.e.,    for  rapid-varying  functions,    the  width     A     should  be 
smaller  than  for  slow-varying  functions. 

Several  graphs  of  the  integrand  of  the  final  equation  (C15)  were 
obtained;  it  is  clear  that  for  a  given  substrate  of  thickness  D  and 
relative  permittivity      Q.     ,    one  can  vary  the  frequency  of  operation, 
width  of  the  strip  conductors  W,    and  separation  between  them  S, 
obtaining  different  corresponding  values  of  effective  wavelengths  on 
the  dielectric     X    .     Therefore,    there  is  no  single  graph  that  could 
be  used  throughout  the  complete  analysis.     However,    a  graph  corres- 
ponding to  a  set  of  parameters  chosen  arbitrarily  will  give  a  fairly 
good  idea  of  how  the  integrand  behaves  and,    consequently,    will 
indicate  the  approximate  limits  of  integration  and  an  acceptable 
segment's  width      A     .     A  typical  graph  is  shown  in  Fig,    5. 

It  was  chosen,    according  to  the  graphical  expression  of  the  inte- 
grand,   that  the  limits  of  integration  will  be  -1700  and  +1700,    since 
contributions  beyond  this  point  were  less  than  1%  of  the  maximum 
contribution.     Likewise,    a  good  choice  for  the  segment's  width 
is  0.5,    which  will  give  a  total  of  6800  steps  of  integration. 
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The  rate  of  change  of  the  value  of  the  integral  in  equation  (C15) 
with  respect  to  the  effective  dielectric  wavelength      A     was  also 
investigated  to  determine  whether  it  was  a  monotonically  increasing 
or  decreasing  function,    an  oscillating  function,    or  a  non-uniform 
function.      For  this  purpose,    the  integral  was  evaluated  at  several 
increasing  values  of     A      and  it  was  observed  that  it  behaved  as  a 
monotonically  decreasing  function,    as  shown  in  Fig.    6. 

The  resulting  integral  in  equation  (CIS), 


i*   ^i  <^(^i  ^  ^''^^  ~(^+w]  sin^  I  In)  dc^ 


■lloo 
will  reach  the  desired  value  of  zero  at  an  approximate  value  of 

for  a  given  set  of  parameters,    i.  e.  ,    frequency,    width  of  conductors, 

separation  between  conductors,    thickness  of  the  substrate,    and 

dielectric  permittivity. 

Since      /o     is  defined  as 


/3 


where     X      varies  as 

it  is  clear  that  an  iteration  procedure  should  be  used,    varying 

from  a  starting  value  equal  to      X      ,    the  dielectric's  wavelength, 

d 

Xd   =T=" 

in  appropriate  steps  toward     A    ,    until  the  value  of  zero  is  reached. 

The  accuracy  of  the  effective  dielectric's  wavelength  depends 
critically  on  the  size  of  the  above  mentioned  iteration  steps. 
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Therefore,    there  is  a  trade-off  in  practice  of  accuracy  vice  computer 
time  or,    equivalently,    accuracy  vice  dollars.      For  the  present  work 
accuracies  of  the  order  of  0.  1%  were  obtained  by  using  iteration 
steps  equal  to  1/1000  of  the  corresponding  free-space  wavelengths. 

The  integration  subroutine  itself  was  tried  for  several  curves 
of  known  equations  of  various  rates  of  variations,    i.  e.  ,    sinusoids, 
paraboloids,    cubics,    etc.     In  all  cases  variations  of,    at  most,    0.2% 
their  theoretical  value  were  obtained.     Also  the  limits  of  integration 
were  increased  to  2500  obtaining  variations  of,    at  most,    0.  05%. 
Finally,    integration  steps  of  0.  4  and  0.  6  were  also  tried  to  check 
accuracy  and,    again,    variations  of  only  0.  1%  were  observed. 

E.      COMPUTER  PROGRAMMING  AND  RESULTS 
1.       Computer  Prograni  Organization 

The  dispersion  characteristic  program,    written  in  FORT- 
RAN IV  language,    accepts  the  following  data: 

-  Dielectric's  permittivity 

-  Width  of  the  conductors,    in  millimeters 

-  Thickness  of  the  substrate,    in  millimeters 

-  Ratio  of  the  thickness  of  the  substrate  to  free- space 
wavelength 

-  Ratio  of  the  separation  between  conductors  to  width  of 
conductors. 
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Likewise,    the  limits  of  integration,    the  step  of  integration        , 
and  the  iteration  step  X,    need  to  be  specified. 

The  program  starts  by  calculating  all  the  required  physical 
constants,    namely,     X     ,    €<,     •     €2    ,    and     jJ©   .     All  the  o«i  -  and 

A -independent  parameters  are  then  evaluated,    i.e.,    k    ,    kp ,    Cd     , 
etc.     The  iteration  step  is  initiated  by  assigning  the  lower  bound 
value  of      A     ,    i.  e.  ,      X    ,    defining  in  this  way  the  value  of    ^   .     At 
this  point  the  integration  is  performed  by  assigning  the  lower  bound 
value  of    o<   ,    i.e.,    -1700  for  the  present  case  defining  consequently 
the  parameters  k^    ,    k^2i    ^c^'    ^   1'  '^  2'    ^^"^    ^   ■^'     ^^ans  are 
provided  to  test  the  real  or  imagimary  character  of    0  2  and  two 
corresponding  paths  then  follow.      The  integration,    or  better,    the 
summation  stops  when    c<    reaches  its  upper  bound  value,    i.  e.  ,    +1700 
for  the  present  case,    and  means  are  provided  to  test  whether  the 

absolute  value  of  the  resulting  quantity  is  less  than  an  arbitrary 

-7 
sinall  epsilon  value,    i.e.,    10        for  the  present  case.     If  the  case  is 

that  the  result  is  greater  than  this  epsilon  value,     A     is  increased 

toward     A     in  one  iteration  step  X  and  the  summation  is  performed 

again.      The  whole  process  ends  either  when  a  result  less  than  epsilon 

.1 

is  found,    in  which  case  the  corresponding      A      is  the  answer  sought, 

or  when      A     reaches  its  upper  bound  value,    i.  e.  ,      A 

The  program  is  outlined  in  detail  in  Appendix  F. 
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Fig.    8. 


Dispersion  characteristic  curves  for  €=16 
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2.       Computer  results 

Effective  dielectric  wavelength  ratio  relative  to  free- space 

wavelength  information  was  obtained  for  three  different  relative 

permittivities,    namely,    for   €.     =   12,    16,    and  20.      This  information  is 

r 

presented  in  terms  of  ratios  of  the  different  parameters  S,    W,    D,    and 
in  order  to  permit  the  use  of  any  set  of  actual  physical  parameters, 
For  each  permittivity,    curves  were  obtained  for  the  following 
values  of  the  ratio  of  separation  between  conductors  to  width  of  con- 
ductors,   S/W  =  %,    %  ,    1,    3/2,    2,    3,    and  4,    corresponding  to  values 
of  the  ratio  of  substrate  thickness  to  free-space  wavelength  D/X 
from  0.  01  to  0.  06. 

Dispersion  curves  for    £      =12  are  presented  in  Fig.    7,    for 

r 

€.      =  16  in  Fig.    8,    and  for    £      -  20  in  Fig.    9. 
r  r 

F.      EXPERIMENTAL  RESULTS 
1.       Fabrication 

The  coplanar  strip  (CPS)  fabrication  was  made  by  photo- 
etching  technique.      Photoresist  is  first  spread  uniformly  across  the 
surface  and  allowed  to  dry.     This  film  is  then  exposed  to  ultraviolet 
light  through  a  photographic  mask  that  permits  the  light  to  fall  only 
on  those  areas  in  which  the  copper  layer  is  preserved.      The  un- 
exposed film  is  next  washed  away  from  the  surface,    leaving  a  solid 
film  of  photoresist  in  the  areas  that  were  exposed.     During  the 
etching  process  the  copper  layer  is   removed  in  the  unprotected  areas 
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but  is  unaffected  under  the  photoresist.      Finally,    the  photoresist  is 
removed  with  acetone,    leaving  only  the  dielectric  surface  with  the  two 
copper  strips,   as  shown  in  Fig.    10. 

It  was  decided  to  construct  the  coplanar  strips  on  a  piece  of 
6"  X  3"  X  1/8"  copper  clad  substrate  of  relative  permittivity  of  12. 
Also,    the  ratio  of  width  of  conductors  to  thickness  of  the  substrate 
was  chosen  arbitrarily  as   1  and,    similarly,    the  ratio  of  separation  of 
conductors  to  width  of  conductors  was  also  chosen  arbitrarily  as   1. 

The  coplanar  strips  were  fed  with  a  transition  consisting  of 
a  2%  inch  piece  of  85  mil  semi-rigid  coaxial  line,  soldered  as  shov/n 
in  Fig.    10. 

2.       Exp er line ntal  data 

Basically,    the  measuring  devices  used  were  a  slotted  line 
fed  by  a  sweep  frequency  oscillator  and  monitored  by  a  voltage  stand- 
ing wave  ratio  meter. 

From  transmission  line  theory,  it  is  known  that  an  open  - 
circuited  line  will  have  voltage  nulls  every  half- wavelength,  which 
in  this  case  will  be  the  effective  dielectric  wavelength. 

Using  this  procedure,    data  was  collected  for  a  range  of 
frequencies  from  1  to  6  GHz  and  a  comparison  was  made  with  the 
corresponding  computer  results,    observing  an  error  of  0.  1%  which 
was,    as  a  matter  of  fact,    the  accuracy  that  was  claimed  for  these 
results. 
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Fig.    10.      Coplanar  strips  fabrication 
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Table  I  lists  the  experimental  results  and  compares  them 
with  the  computer  results. 
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TABLE  I 


DISPERSION  CHARACTERISTIC  COMPARATIVE  RESULTS 


Coplanar  strip  of  characteristics   :^^  =  12 

d  :=  0.  125" 
s  =  w    =  d 


Frequen 

icy 

^/X     from 

^/\     from 

GHz 

Fig.    7 

experiment 

1.0 

0.447 

0.446 

2.0 

0.435 

0.435 

3.0 

0.425 

0.424 

4.0 

0.415 

0.414 

5.0 

0.406 

0.406 

6.0 

0.398 

0.399 
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III.  CHARACTERISTIC  IMPEDANCE 

G.   THE  AVERAGE  POWER  AS  A  FUNCTION  OF  THE  HERTZIAN 
VECTOR  POTENTIAL  FUNCTIONS 

As  stated  by  Collin  [Ref.    4],    a  general  expression  for  the  tirae- 

average  power  flow  in  terms  of  the  electric  and  magnetic  fields  is 

"Rv.  ('^.y)  =  i(?e [lis '<'?*•  aj  ds]  <Gi) 

where, 

(G2) 


%x>W*.a2  -  <i^^y*-<§yvyx    . 


Referring  back  to  equations  (Al)  through  (A7),    it  is  found  that  the 
equations  for  the  electric  and  magnetic  fields   can  be  expressed  in  terms 
of  Hertzian  potential  functions  as, 


a        /v  c)S>       •  dot 


(G4) 
(G5) 
(G6) 


1     ^  (G7) 


^r(^^--3-^) 


<^       "  ^^    '  (G8) 
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So,    substituting  terms  in  the  expression  for  average  power, 
^  (".Yl  4  "^e  [  ji  (  ix  Jiy*  -  fey  '"''* )  M 

Taking  the     OC  -domain  Fourier  transform  of  the  average  power 
expression, 

7  /  (GIO) 

where  the  Fourier  transform  of      ^/6\  has  been  applied,    i.  e.  , 

Frora  equation  (GIO),    performing  the  operations  indicated, 

"^-EKy)  =  ^ ^  [  If  l*"'/^^^  ^'^M^  ^''(«.V) -^  '^^/^  ^''^^^'Vi ^'"fM^ 


(Gil) 


As  far  as  the  first  four  terms  in  the  above  expression  for 
average  power  are  concerned,    they  will  be  real.     However,    for  the 
remaining  four  terms,    since  they  correspond  to  cross-power  terms, 
a  more  careful  analysis  has  to  be  done.     A  detailed  exposure  to  such 
analysis  is  given  in  Appendix  G. 

The  real  part  of  the  average  power  turns  out  to  be. 
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+•# 


5-  (-'i)  =  2^  \\  {-^/^  (^^- V)(e-  An.^  ^/J.|f.^V.)n  e 


2.   -2r,(Y-<<) 


2^ 


^a)^(^^->iO(€.  X)^^V^^  ^/j.  \'b''^(^)\ )  e^''1 


(G13; 


-f  COS 


+  ^Sc,'  Vj  (iB^V-lIc'^M  -  B'^M  |c«m|)]  I 

I 


and  similarly,    for    Y       imaginary, 


4« 


,.u>/3  (^^>r/)(€.  ^^4^.  |x)^'^cmI')  e'''''/ 


-c< 


(G13a) 
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It  is  clear  from  the  above  expressions  that  one  can  integrate  out 

the  y-dependence  of  the  average  power,    leaving  only  the  oC  -dependence. 

The  integration  details  are  shown  in  Appendix  G,    the  resulting 

expressions  being, 
4«e 


^AOt 


+  •1 


and,    similarly,    for     Y       imaginary, 


"^t  (-1  = 


(G14) 


(G14a) 
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Fig.    11.     Z  -  direction  surface  current  density 
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4-p.(lf^'\.)\'^\'D'^^U)\')] 


where  in  both  equations  (G14)  and  (G14a),    the  equality    Y     =  V     , 

1  3 

has  been  applied. 

H.      THE  CHARACTERISTIC  IMPEDANCE  IN  TERMS  OF  THE 
DISPERSION  CHARACTERISTICS 

In  general,    one  could  express  the  average  power  in  terms  of  the 

current  in  the  coplanar  strips  as, 

9^,^il'^'  (HI) 

from  where, 

p   ^2:^  (H2) 

In  the  present  work,    it  was  assumed  in  the  dispersion  character. 
istic  part  that  the  surface  currents  through  the  strips  were  uniform; 
furthermore,    the  current  can  be  expressed,    as  shown  in  Fig.    11,    as 

strlj, 
where  the  approximation  is  valid  since  uniform  surface  current 

densities  were  assumed. 


46 


Furthermore,    since  it  was  also  assumed  a  unit  amplitude  square 

pulse  as  the  axial  surface  current  density,    i.  e.  ,    J    (x)  =  1,    it  follows 

z 

that, 

(H4) 
I  =  W  I 

I^=  W^    .  (H4a) 

All  constants  in  the  average  power  expressions  of  equations 
(G14)  and  (G14a)  are  completely  defined  by  choosing  appropriate 
values  for  the  width  of  the  strips,    separation  between  strips,    thick- 
ness of  the  slab,    dielectric's  permittivity,    and  frequency  of  operation. 

The  propagation  constant    A  is  now  well-defined  from  the  value 
of  the  corresponding  effective  dielectric's  wavelength     X    ,    i.e., 

'         X'  (H5) 

where    X    is  determined  from  the  dispersion  characteristic  curves  of 
Figs.    7,    8,    and  9. 

Therefore,    the  expressions  to  be  evaluated  for  CPS  characteristic 
impedance  are, 


(H6) 


and,    similarly,    for  the  case  where     o         is  imaginary, 
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I.       NUMERICAL  INTEGRATION 

The  numerical  integration  problem,    as  far  as  the  characteristic 
impedance  is  concerned,    is  the  evaluation  of  equations   (H6)  and  (H6a). 
Again,    although  the  limits  of  integration  in   o^  are  -  oo  and  +  oo,    a 
prior  knowledge  of  the  behaviour  of  the  integrands  allows  the  use  of 
truncated  limits. 

As  in  the  dispersion  characteristic  part,    the  numerical  integra- 
tion method  used  is  a  modified  Simpson's  rule.      Practical  truncation 
limits  are  again  -1700  and  +1700  in  steps  of  c<  =  0.  5,    which  will  give 
6800  computing  points  for  a  given  set  of  parameters. 

Once  the  physical  parameters  are  specified,    namely,    the  sub- 
strate's relative  permittivity,    the  separation  between  conductors  to 
width  of  the  conductors  ratio,    and  the  substrate's  thickness  to  free- 
space  wavelength  ratio,    all  the  constants  are  well-defined  except  for 
the  propagation  constant    A  .     However,    this  value  can  be  interpolated 
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from  the  dispersion  characteristic  curves  shown  in  part  E.      There- 
fore,   numerically- wise,    the  calculation  of  the  characteristic  imped- 
ance is  simpler  than  that  for  the  dispersion  characteristics  since  no 
iteration  procedure  is  required. 

J.       COMPUTER  PROGRAMMING  AND  RESULTS 
1.       Computer  program  organization 

The  characteristic  impedance  program,    written  in  FORT- 
RAN IV  language,    accepts  the  following  data: 
Dielectric's  relative  perinittivity 
Width  of  the  conductors,    in  millimeters 
Thickness  of  the  substrate,    in  millimeters 

Ratio  of  the  thickness  of  the  substrate  to  free-space 
wavelength 

Ratio  of  separation  between  conductors  to  Vv^idth  of 
conductors 

Ratio  of  effective  dielectric's  wavelength  to  free-space 
wavelength 

Likewise,    the  limits  of  integration,    and  the  step  of  integra- 
tion need  to  be  specified. 

The  program  starts  by  calculating  all  constants  and  para- 
meters and  directly  evaluates  the  impedance  integral;  basically,    the 
procedure  is  the  same  as  in  the  dispersion  characteristic  part  except 
for  the  iteration  search  for     X    since  now  it  is  completely  defined. 

The  program  is  outlined  in  detail  in  Appendix  H. 
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2.       Comp-uter  results 

Characteristic  impedance  information  was  obtained  for  three 
different  relative  permittivities,    namely,    for    G  „  =  12,    16,    and  20. 

As  in  the  dispersion  characteristic  part,    this  information  is 
presented  in  terms  of  ratios  of  the  different  parameters,    i.  e.  ,    S/W 
and  D/  X     ,    in  order  to  permit  the  use  of  any  set  of  actual  physical 
parameters. 

For  each  pernnittivity,  curves  were  obtained  for  values  of 
S/W  equal  to  %,  %,  1,  3/2,  2,  3,  and  4,  corresponding  to  values  of 
D/  X      from  0.  01  to  0.  06. 

Characteristic  impedance  curves  for    ^      =12  are  presented 
in  Fig.    12,    for    e       =   16  in  Fig.    13,    and  for    e      =  20  in  Fig.    14. 
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Fig.    12.      Characteristic  impedance  curves  for  &=  12 
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IV.     APPLICATIONS 


The  present  method  is  applicable  to  problems  based  on  any 
coplanar  strip-type  configxiration  and  also  to  problems  related  to 
slot  lines,    such  as  coupled  slots  or  coplanar  waveguides. 

For  the  latter  case,    the  configuration  is,    in  general,    as   shown 
in  Fig.    15,    and  the  main  difference  with  the  coplanar  strip  problem 
is  the  knowledge  of  the  electric  fields  instead  of  the  surface  current 
densities. 

Referring  back  to  equations  {B46)  and  (B47),    these  can  be 


expressed  as, 


3iM 


%%(ti)      M(»(fc<,(6) 


where, 

therefore, 

Similarly, 

M,  (ti.^  a)        £?  M 


Jk('<^ 


(Li: 


(L2) 


(L3) 


(L4) 
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and, 

In  general,    for  any  slot-type  configuration,    z-directed  electric 
fields  can  be  neglected  and  the  transverse,    or  x-directed  electric 
fields  can  be  approximated  with  square  pulse-type  distributions. 

Unlike  the  coplanar- strip  problem,    here  it  is  possible  to  obtain 
an  odd  and  an  even  type  of  electric  field  excitations  across  the  slots. 

In  either  case,    equations  {L3)  and  {L5)  reduce  to,    respectively. 


*^U\f>) 


^xi(^)  =  J^^*^^  (L6) 


and 


■;;f7^§K,(^c^-^x(.)  (1.7) 

where  i  =  e,    o  depending  on  the  type  of  excitation,    even  or  odd,    as 
shown  in  Fig.    16. 

As  far  as  the  odd  case  of  excitation  is  concerned,   ^j^(<*:l assumes 
a  form  similar  to  J   loc)  in  equation  (Gil),    i.e., 

<§Xo  ("<)  =  ^  ^'^  T  (s+vo^  s/o  I  tvj  (L8) 

°  tX.  2.  ^  2.  . 

Performing  the  necessary  Fourier  transform  and  algebraic 
reductions  becomes,    for  the  even  case  of  excitation, 

%>Kc  (^)  =  ^  ^  I  f^-^^^  ^'<r^f\^},  (L9) 
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Applying  Galerkin's  method  and  by  Parseval's  theorem,    in  the 
same  approach  as  for  the  coplanar  strip  problem,    the  final  iteration 
equation  for  the  effective  dielectric's  wavelength  becomes,    considering 
only  equation  (L7), 

J     "^JZJ^t^    ^Xc(^)  fexo^^^   cJt^^CJ  (LIO) 

or, 

where  M   {^  ,  f?>)  through  M    (  «<  ,  ft  )  are  defined  in  equations  (122) 
through  (125),    respectively. 

Remaining  procedure  is  correspondingly  similar  to  that  for  the 
coplanar  strip  case. 
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Fig.    15.      Coupled  slots  configuration 
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Fig.    16.     Even  and  odd  transverse  electric  fields 
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V.     CONCLUSIONS 

The  method  used  in  this  present  work  succeeded  in  providing  an 
accurate  set  of  results  for  the  dispersion  characteristic  and  the 
characteristic  impedance  of  a  Coplanar  Strip  transmission  line. 

It  was  also  discussed  that  an  extension  of  this  method  will  allow 
corresponding  results  for  the  Coupled  Slots  or  Coplanar  Waveguide 
configuration.     However,    for  the  latter  case,    care  must  be  taken  in 
adapting  the  equations  to  the  even  case  of  excitation  due  to  the  fact 
that  the  Fourier  transform  of  the  transverse  electric  field  will  now 
be  real  in  character,    instead  of  imaginary  as  in  the  odd  case  of 
excitation. 

The  amount  of  algebraic  work  can  be  reduced  substantially  by 
the  use  of  matrices,    which  for  the  present  configuration  will  have 
dimension  eight.     However,    the  algebraic  approach  used  here  gives 
much  more  insight  to  the  problem  and  also  reduces  sharply  the 
required  computer  time. 

As  far  as  the  theoretical  development  is  concerned,    the    only 
approximation  used  in  the  representation  of  the  axial  surface  current 
density,    i.  e.  ,    a  square  pulse,    and  the  value  of  zero  assumed  for  the 
transverse  surface  current  density.     Depending  on  the  desired  accu- 
racy,   more  complete,    and  indeed  complex,    representation  of  these 
surface  current  densities  can  be  used. 
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In  the  numerical  integration  part  of  the  problem,    the  integrating 
method  itself  has  been  proven  to  be  sufficiently  accurate.      In  fact, 
the  only  significant  approximations  used  are  the  choice  of  the  trun- 
cation limits,    the  integration  step  and  the  iteration  step.     However, 
several  different  sets  of  these  choices  were  tried  and,    at  least,    for 
the  present  configuration,    negligible  variations  were  observed. 

It  is  clear  that  the  two  programs  used,    i.  e.  ,    the  dispersion 
characteristic  and  the  characteristic  impedance,    could  be  lumped 
together  into  one  program  reducing  in  about  half  the  araount  of  compu. 
tation  and,    therefore,    the  required  computer  time. 
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APPENDIX  A 
AUXILIARY  VECTOR  POTENTIAL  FUNCTIONS 

A.      MAGNETIC  HERTZIAN  VECTOR  POTENTIAL  FUNCTION 

Consider  a  homogeneous,    source-free,    isotropic  region;  hence, 
there  is  no  charge  density  and 

V-S  =  o  .  (1) 

From  Vector  Algebra,    it  is  known  that  the  divergence  of  a  vector 
is  zero  if  the  vector  is,    in  turn,    the  curl  of  another  vector;  there- 
fore,   one  can  state, 

:  ^-^^^i,  =o       >^Tt\,  (2) 

^--jco^VxTCi,      ^  (3) 

For  time -varying  fields 


so  applying  equation  (3), 


(4) 


(6) 


where. 

Similarly,    for  time-varying  fields, 
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so,    applying  equation  (3), 
Vx(~jc0^^x7(^)  =-^co^(^%  i-t^CJ))  (8) 

(9) 
Up  to  this  point  0  has  not  been  defined;  arbitrarily  choose 


(f)  =  \7'K^ 


(10) 


Then,    equation  (9)  becomes, 
Also,    it  is  known  that, 

Summarizing,    for  TM  modes, 

Similarly,    for  TE  modes, 

>^.  oi.   ^- 
«  v7k\7x)I^ 
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B.   TE  AND  TM  MODES  FROM  VECTOR  POTENTIALS 
1.   TE  Modes 


=1  Vx'^xTL*' 


For  TE  modes,       ^        =0;  therefore, 
from  where, 


It-^-^a,.^ 


(17) 
(18) 


where     It**     satisfies: 
V2"ft,H^Jl^Tt»'  »0.  (19) 


(20) 


TL^-rJli-O,  (21) 

Let: 

K^=  <;6^c-^^a^  .  (22) 

Then,    froin  equation  (19)i 
Vt(p^'^^4>^=^o  (23) 


(24) 


where. 

Summarizing, 

(27) 
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2.       TM  Modes 

For  TM  Modes,    K        =0;  therefore, 

z 


(28) 
(29) 


(30; 


and, 

^t^<^%t?0^=O  (31) 

where, 

Summarizing, 

^  (33) 

(34) 
^""X  "Y'  ^^  ^  (35) 

Equations  (Al)  and  (A2)  correspond  to  equations  (33)  and 
(25),    respectively. 
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APPENDIX  B 
TRANSVERSE  ELECTRIC  AND  MAGNETIC  FIELDS 

Consider  Maxwell's  point  form  of  Ampere's  law, 

"$-■;;  «ja3€l  (36) 


from  where, 


7 


a-e 


-f  =  j"^e  ^x 


apply, 

obtaining  from  equations  (38)  through  (40), 


d 


y^K  -    TT    ^J**^^   ^ 


'^-  f  -  J-  ^^  . 


(37) 


(38) 


a^y       Mx^.,...   ;£.  (40) 


and,    dvxe  to  the  z-dependence  of  the  fields,    i.  e.  ,     C        ,    one  can 


(41) 


-y^y    -  ja>e  ^x  (42; 

a>x     o         7  (43) 


(44) 
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Similarly,    starting  frora  Maxwell's  point  form  of  Lenz'  law, 
from  where, 

Substituting  the  value  of     (o      from  equation  (48)  into  equation 


X 


(42),    one  obtains, 


f-v«/=f  (t-J-/>^r) 


1 

^Hi     .    ,      d$i  (50: 


Similarly,    substituting  the  values  of    <S       from  equation  (47) 
into  equation  (43),    of  9^      from  equation  (43)  into  equation  (47),    and 
of    9^       from  equation  (42)  into  equation  (48),    one  obtains, 

y 

-   For   ^     , 

X 
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SUi      .     .    f   /  ^^     . 


7 


«.(vw^.)-^lr-j-'^  (51) 


-  For  %     , 

-  For    %     , 

X 

Summarizing, 


(54) 
aee^. a^a\  (55) 


Equations   (A4)  through  (A7)  correspond  to  equations   (54) 
through  (57),    respectively. 
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APPENDIX  C 
DERIVATION  OF  FIELD  EQUATIONS 

The  differential  equation  defining  the  field  functions  is,    as 


expressed  in  equation  (B5), 

In  turn,    the  propagation  constant  X ■      is  expressed  as. 


-^2   ^i(^,<j)-yL'^i('^.'^)^0     ,       6.1,-2,3.  (58) 


As  stated  before,    the  propagation  constant  for  air  medium, 
regions   1  and  3,    is  always  real.     However,    for  the  dielectric  medium, 
the  propagation  constant  will  be  imaginary  for  small  values  of   <^   and 
real  for  large  values  of  o<  . 

For  region  1,    the  field  function  differential  equation  solution  is 
as  follows, 

Clearly,    for  infinite  positive  value  of  y,    the  field  should  vanish, 
or. 


(6i: 
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Similarly,    for  region  3, 

In  this  case,    for  infinite  negative  values  of  y,    the  field  should 
vanish,    or, 

f-"  ^  (63) 

Finally,    for  region  2,    there  are  two  solutions,    corresponding  to 
the  real  or  imaginary  character  of  the  propagation  constant;  the  solu- 
tion corresponding  to  the  imaginary  case  will  first  be  sought  as 
follows, 

For  values  of    c<,  greater  than    ^  V^r-J    ,    the  propagation  constant 
will  be  real,    or 

i>^(b<,ij)  ^^fS.in\^  Y^u '¥)C;tCx:>6,\n  Yi^Q  (65) 

where, 
^T'-j"^  (66) 

(67) 
One  can  define  the  following  constants, 
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AH  =  IC2 

B,  (u)  ==  Ky 

then,    finally,    applying  the  functions  to  the  electric  and  magnetic    field 

cases, 

5/'V^,Y)=A^"(^)  «"'''''  (69) 

^^'('•c,*/)  «  S^*'^'^)  si^A  V^^  +  C^^Vv)  co5li  Ve^Y  (70) 

^^\c<,^)  =  bJ^Vm  si/?  Y^tJ  ^  C^*Vk)  cdS  n,y  (70a) 

^fc^.^)^  I>''hu)e''^1  (71) 


Equations  (69)  through  (74)  correspond  to  equations  (BIO)  through 
(B15),    respectively. 


(72) 
(73) 
(73a) 
(74) 
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APPENDIX  D 
MATRIX  CURRENT  AND  FIELD  EQUATIONS  DERIVATION 

From  equations  (B32)  and  (B33),    one  could  express  the  constants 
A^      '(K)andA         (  C<  )  in  terms  of  C^      '(  (x  )  and  C*      '(«)as, 

ft<"(«,  -  {^f{ c"'m [^^  (^y  s/oJ,  y. d  ^  cosh  y,  j] 


a'^'(.>  = 


Define  the  following  constants: 


Therefore,    equations  (75)  and  (76)  become, 

A"'  (k^  =  o,  c^^'ck)  ~  Oj.  a}^^  (t,) 


(77) 
(78) 
(79) 
(80) 


(81) 
(82) 
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from  where  one  could  solve  for  c'       '{^)  and  C^        (  C<  )  in  terms  of 
a(  ^  )(oc)  and  A^      ^  (X  )  as  follows, 


(83) 
^^\^.  -  i  r.v    A<*V.x     .,    Aft^^.n  (84) 


Once  again  define  the  following  constants 


R= ^^i^-—  (85) 


^, !Ji (86) 

K)|fJn~N)2.K)A 

^-  *^3>  (87) 

P  n3/  (88) 

Therefore,    equations  (83)  and  (84)   can  be  expressed  as, 


which  correspond  to  equations  (B42)  and  (B43),    respectively. 

Substituting  equations  (89)  and  (90)  into  equations  (B34)  and  (B35), 
one  obtains, 


(91) 


-t^.^-Vc^fC-   [(^)^-j]^.|(^i)^P^).;,,,,. 
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--i-^  ^  (-fe)'  P.  ]  -  --t.  n ^  [J-.r,(^: )' P,  ,92) 

Define  the  following  constants: 

(93) 


-  (oeoY, 


(96) 
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Therefore,    equations  (91)  and  (92)  beconae, 

(  ^  )  (  ^^  ) 

from  where  one  could  obtain  expressions  for  A         {  ^  )  and  A         (  ck  ) 

in  terms  of  J   (^  )  and  J    ( c»<  )  as  follows, 
X  z 

'^'*'^"^^;; — ^7r;7-lG?.Je(^i-Q3JvC^)]  (99) 

A*'^M  =  - L__[Q^J^(^)+6?W,(^^)]  (100) 

which  correspond  to  equations  (B44)  and  (B45),    respectively. 

Substituting  equations   (99)  and  (100)  into  equations   (B36)  and  (B37), 
and  solving  for     ^     (  oc  )  and    <o   (      )  in  terms  of  J   {^  )  and  J   {c>c  ),    oj 


)ne 


X 

obtains: 


^ <2.  J^^O-QsJxC*^)      =    4iM 

(-J^yIq.  Ji(M-G?3Jx(«)l 

C?,Q4+6?zQa    I  ^        1.    '     "^  J  (102) 

Define  the  following  constants: 

M,  =  - 

QtQ^-i-QzQs,  (103) 


Mz 


-'tc,'^<Qa 


G?,G?^ -t-Qz^Ja  (104) 


(106: 
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Then,    finally,    equations   (101)  and  (102)  become, 
M,KfS^  J^(v)   4    Mz(k,|S)  OxC*^)    =    ^;^M  (107) 

M3(^/2>)  J£(»<)  +  M4  (o<,/2>)  0x(6<)   ^    %x(^]  (108) 

which  correspond  to  equations   (B46)  and  (B47),    respectively. 
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APPENDIX  E 
CURRENT  AND  FIELD  EQUATIONS  RE-DEFINITION 

From  equation  (A3),    the  character  of  the  propagation  constant 
was  defined  as  lossless,    i.  e.  , 

Y-j(5   .  (109) 

So,    equations  (77)  through  (80)  are  modified  as, 

-^^  (fen  Site) -•^^'''^-^^^<^i  (no, 

'^-j(fe]VSrKfe)^-0-''V-.j  (111) 


Consequently,    equations  (85)  through  (88)  become. 


K).>J4-*-ll04K)al 

IkJ^I 

kJ.kJu+IOzJkJsI 

Ns 

fOiKJa  +  lOzIkIa 

k3, 

(113) 


(114) 
(115) 


^^"J  .>  .^■'.7!  If. 77  (116) 


(117) 
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(118) 


(119) 


Similarly,    equations  (93)  through  (96)  are  re-defined  as, 

Finally,    equations   (103)  through  (106)  become, 
M,-tl — i i — r  (122) 

^'j     lQzlQ3-Qil«4l  (123) 

M^j  ; i r  (^P^'^i  iQzl-^^fS^S') 


(124) 
(125) 


Furthermore,    since  from  equation  (Cll)  J    (  ^  )  is  imaginary  and 

z 

J    (^  )  is  zero,    the  expressions  for  the  different  constants  A  through 

X 

D  need  to  be  re-defined. 
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Equations  (B44)  and  (B45)  become, 

which  cori*espond  to  equations  (C16)  and  (C17),    respectively. 
Similarly,    equations  (B42)  and  (B43)  become, 

C'W  (^1  =j  [|P4  a'"m  -?,  |A*V«)|]  (129) 

which  correspond,    in  t-urn,    to  equations  (C18)  and  (C19),    respectively. 
Finally,    equations  (B38)  through  (B41)  become, 

^"'m  -  (^;f  (L'Vi  (132) 

which  correspond  to  equations  (C20)  through  (C23),    respectively. 
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APPENDIX    F    -    DISPERSION    CHARACTERISTIC    PROGRAM 

THIS    PROGRAV,    SOLVES    THE    DISPERSION    CHARACTERISTICS    OF 
COPLANAR.  PARALLEL    STRIPS    ON    A    DIELECTRIC     SUBSTRATE. 

PROGRAM    ASSUMES    FOLLOWING    PARAMETERS    ARE    K^CWN 

-  THICKNESS    OF     THE     SLAB, DD, IN    MILLIMETERS 

-  WIDTH    OF    THE    STRIPS, WW,  IN    MILLIMETERS 

-  DIELECTRIC'S    RELATIVE     PER  MI TT I VI TY , E R 

-  RATIO    OF    SEPARATION    BETWEEN    CONDUCTORS    TO    WIDTH    OF 
CONDUCTORS, SOW 

-  RATIO    OF     SUBSTRATE'S    THICKNESS    TO    FREE-SPACE 
WAVELENGTH,  DRAT  10 

THE    METHOD    USED     IN    FIMDING    THE     EFFECTIVE     WAVELENGTH    IN 
THF    STRUCTURE, LPRIME,  IS    AN     ITERATION    METHOD, WHICH 
SOLVES     FOR    THE    ZERO    VALUE    OF    AN     INTEGRAL     IN    THE    ALFA- 
DOMAIN    FOR    PROGRESSIVELY    LARGER     VALUES    OF    EFFECTIVE 
WAVELENGTH. 

THE    SUMMATION     IN    ALFA    IS    AN    APPROXIMATION    TO    AN 
INTEGRATION    IN    THE    ALFA-OOMAIN    Fi^OM    (-INFINITY)    TO 
(+INFINITY). 

THE    APPROXIMATIONS    USED    IN    THE    PRESENT    PROGRAM    ARE 

-  THE     ITERATION     STFP,X 

-  THE    EPSJLON    VALUE, DELTA 

-  THE     INTEGRATION    STEPrA 

-  THE    LIMITS    OF     THE    ALFA-DOMAIN    I NTEG R AT  I  ON , ALFA 
AND    B. 

IT    IS    CLEAR    THAT    ANY    OF    THESE    VALUES    MAY     EE    CHANGED 
ACCORDING    TO    THE    SPECIFIC    CASE    GR    DESIRED    ACCURACY. 

THE    QUANTITY  EO    IS    THE    FKEE-SPACE    PERMITTIVITY    IN    MKS 

UNITS. 

THE    QUANTITY  MU     IS    THE    FREE-SPACE    PERMEABILITY    IN    MKS 

UNITS. 

THE    QUANTITY  C     IS    THE    FREE-SPACE     SPEED    OF    LIGHT    IN 

METERS    PER    SECOND. 

THE    QUANTITY  Ml     IS    THE    STRIP    CURRENT    IN    AMPERES. 

THE    VECTOR    SOW    STORES    THE    RATIOS    OF    SEPARATION    OF 
CONDUCTORS    TO    WIDTH    OF    CUNOUCTOkS    WHICH     APE    DESIRED 
TO    SOLVE    FOR.    THE    NUMBER    J    SETS    THE    DIMENSION    OF     THIS 
VECTOR. 

PROGRAM    DEVELOPED    BY    LT(JGi     ARMANDO    LUNA    ECHEANDIA, 

PERUVIAN    NAVY. 

SUPERVISOR    -     PROF.     JEFFREY    D.    KNORR,PH.D. 

U.S.     NAVAL    POSTGRADUATE    SCHOOL, MONTEREY, CAL IFORNIA 

SEPTEMBER, 1973. 

IMPLICIT    REAL*^4     (M,N,K,L) 
DIMENSION    SOW (7) 

DATA    WW/3.17/,DD/3.17/,ER/12.0/ 
DATA    A/0.5/ 
DATA    B/1700./ 

DATA    X/.4E-4/,0ELTA/1.0E-7/ 
J=7 

PI=3. 14159 
E0  =  l.E-9/(36.'S^PI  ) 
MU  =  4.E-7-?PI 
C=3.E8 
FACT0R=1.E9 
W=WW*1. E-3 
D  =  DD''l.E-3 

READ  ( 5,100)  (SOW( I )  ,1=1  ,  J) 
100  FORMAT  (7F5.3) 
24  READ  (5  ,18)  DRAT  10 
18  FORMAT  (F6.4) 
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c 
c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
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FPEQ  =  DRATIO-C/D 
FRE=FREQ/FACTOR 
LPRIME=C/(FREQ-SQRT{ER)  ) 
0NEGA=2  ."^PI     FREQ 
0E  =  OMEGAi<E0'  ER 
OM  =  niMEGA-^-MU 
K1=CMEGA^-5QRT(MU^E0  ) 
K2=K1-SQRT(ER} 
K3  =  K1 
Kll=Kl--^2 
K22  =  K2^  ^-2 
DO    101     I=1»J 
S=SOW( I ) 
22    BETA  =  2.-'PI/LPRIME 
B2=BETA^"2 
KC1SQ=K11-B2 
KC2SQ=K22-B2 
KC3SQ=KC1SQ 
KC4=(KC2SQ/KC3SQ) 
KC5=KC4-1.0 
KC6=KC5^'^  2 
M3=0.0 
ALFA=-17C0. 
1     ALPhA=ALFA+^/2.0 
A2  =  ALPHA-^t.2 

PHYSICAL    CONFIGURATION    DEPENDENCE    AS    DERIVED    BY    THE 
FOURIER    TRANSFORM    OF    THE    CURRENT    DENSITIES. 

F1  =  ALPHA-^  (S  +  W)/2  .0 
F2^ALPHAv W/2.0 
S1=SIN(F1) 
S2=SIN{ F2) 

53  =  Slv:^*2 

54  =  S2=«^^-2 
AB=ALPHA' BETA 
G1=SQRT{A2+B2-K11) 
G3  =  G1 
GI1=GI^''2 
G22=A2+B2-K22 

IF    (G22.LT.0.0)     GO    TO    7 

HYPERBOLIC    CASE,     GAMMA    2    IS    REAL 

G2=SQRT{G22) 
G2D=G2-D 
SI=SIfSlH(G20) 
C0=C0SH(G2D). 

INTERMEDIATE    ALGEBRAIC    STEPS 

N1  =  KC4^  (G3/(G2'^£R  )  VKC4-SI+C0) 

N2  =  KC4'  AB/(OE-'-G2)^Kr5'-^SI 

N3  =  AB/(0M>G1)>KC4^  G3 /  (  G2-'- ER  )  ^  KC  5^  SI 

N4  =  l./  (0M^^G1)>^  (  (OM'  G2-(ABV-^2)--KC6/(0E-G2)  )"  SI  + 
10M^G3-^KC4+C0) 

N5=N1=^N4+N2*N3 

P1=N4/N5 

P2=N2/N5 

P3=N3/N5 

P4=N1/N5 

Q1  =  KC  2  Sg>'  (  (  G3  /G2  )  -  KC4-'  P4-AB/  (  CM-1G2  )v  P2  •■  KC5-'  S  I  +KC2  SQv 
IPA^CO+KCISQ 

Q2=KC2S0-  (  (AB/(0M->G2}-'KC5^  P1+(G3/G2)=  KC4--P3)^  SI+P31C0) 

Q 3=  (  AB  -  (  G3  /G 2  )  ^  KC  4<  P4-0E  ^  G2»  P2-  (  AB-t  '  2  )  /  (  CK-  G2  )  =•  KC 5  '^ 
lP2)*SI-(  (OE/FP  )--G3'KC4--P2-AB-P4-AB<P4'KC5)^^Ca  +  AB 

Q4  =  ( (  AB  ^:-  *  2  )  /  (  CM^  G2  ) '  KC5  -  PI  +  A  B  -  (  G3  /  G2  )  -'  P3  :  KC4+0  E  ^'  G2-' 
1P1)*^SI+(A3  ;P2-K0E/ER  )'-G3^KC4    P1  +  A6  ^-KCS^^  P3  )  ^v  C0+ (  OE/ER)-^ 
2G1 

GO    TO    9 

TRIGONOMETRIC    CASE,     GAMMA    2     IS     IMAGINARY 
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c 
c 
c 


G22=ABS (A2+B2-K22) 

G2=SQRT(G22) 

G2D=G2^D 

TSI=SIN(G2D) 

TC0=C0S(G2D) 

INTERMEDIATE    ALGABRAIC    STEPS 


Nl  =  KC4r 
N2=^KC4--^ 
N3=AB/( 
N4=l./( 

lOM^  G3^^K 
N5  =  N1^N 
P1=N4/N 
P2=N2/N 
P3=N3/N 
P4=NI/N 
01=KC2S 

1KC2SQ-^P 
Q2=KC2S 

1P3^TC0) 
Q3=(AB-'^ 

1P2}'^TSI 
Q4=(( AB 

IPD-^TSI 

2(CE/ER) 
9  Q5=Q2^^^Q 
M1=KC1S 
M2  =  16.' 
M3=M3+M 
ALFA=AL 
IF  (ALF 
GO    TO    1 


(G3 
AB/ 

GM^ 

4+N 
5 
5 
5 
5 

0-  { 
4vT 
0-.  ( 


/(G2'  ER)^'  KC4-TSI  +  TC0) 

(OE:'  G2}^KC5  'TSI 

G1)^KC4'<G3/(G2-!  ER)  '  KC5=  TSI 

Gl )  =>••  (  (  -OM f  G  2-  ( A  B * ^  2  )  >•■-  KG  6/  ( OE  f  G 2 )  )  -  TS I  + 

TCD) 

2vN3 


(  G3  /G2  )  -KC4--  P4- AB/  ( OM^  G2  J*  P  2  ■  K  C  5  ) '  T  S  I  + 

CC+KCISQ 

(Ab/(GM--G2)  -KC5  ■P1  +  (G3/G2}  .KC4rp3  )rTSI  + 


(G3 
-(  ( 

+  {  A 

^  Gl 

3-Q 

G*Q 

Ml./ 

2-' A 

FA  +  A 

A.GT.B  )    GO    TO    2 


/G2  )  -^KC 4 -  P4+G E^- G2-^^  P2-  (  A B^'.'2  )  /  (  ai/-^^G2  ) ■■  KC  5^= 
OE/  ER)-r;G3KC4'  P2-AB-P4-AB  ■P4:-KC5)- TCC+-AB 
)/(nM'-G2)^  KC5    PH-AB^-  (  G3/G2)-  P3 -  KC4-0E    G2* 
B-'P3+(  OE/ER  )-G3^KC4^  P  1+AB    KC5^  P3)*  TCO+ 

1^  Q4 
1/Q5 
'A2VS34  34 


ITERATION    STEP 

2     IF     (M3. LT. DELTA) 
LPRIME^LPRIME+X 
GO    TO    22 
19    WRITE     (6, 
191    FORMAT( • 

1»AND    FREQ 
101    CONTINUE 
GO    TO    24 
23    STOP 
END 


191  ) 
'  -/, 


GO    TO     19 


S,FREQ, LPRIME 
3Xt'LAH3DA    PRIME 
»E16.7,2X,MS     S 


FOR    S    = 
E  1  6 .  7 1  / ) 


• ,F10.5»2X 
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APPENDIX  G 
CHARACTERISTIC  IMPEDANCE  RELATED  PARAMETERS 

It  is  convenient  to  recall  the  Hertzian  potential  functions  for  the 
three  regions  of  analysis.      The  electric  field  functions  are, 

i,,«i(.,Y)=A«'Me-'^'V-*"  (134) 

*f' K^l^  8'%sm(i  r^u  +  <!'"c«<)  toiU  Yjy  (135) 

C<«<.f)'D'Ve'''y  (136) 

from  where  the  corresponding  y-derivatives   can  be  obtained  as, 
-'     =-Y,A^^Vo<\€        '  (137) 


^ 


1 


,<c) 


a^r   ..  ^^cc)..      ....  .«r) 


—  2 


=  Y2  [b^''V»<)  tosh  V;.<^  +  C^'V^)6;^Ky;,tv]  (138) 


'1 

^"=Y.I>%e'^^  (139) 

As  stated  in  the. dispersion  characteristic  part  of  the  present 
work,    specifically  in  equations  (C16)  through  (C23),    the  magnetic 
field  related  constants  have  imaginary  character.      Therefore,    it  is 
convenient  to  modify  the  magnetic  field  Hertzian  potential  functions 
to  account  for  this  fact. 

3&f  c«,.f)»j  [(b'"(«)|  iif<ur,u  +  (c*V.»|Qjsh  nw]  ,j4i) 
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and,    again,    the  corresponding  y-derivatives  can  be  obtained  as, 

3_=-jr,  lA^^^Klje       ^  .  (143) 

^  =j'>'^[|B^^V^)l  c^K  V,c^+|c<^VK)lsioh  Yzcj]  (144) 

^jL^jY^ll^^^l/^y  (145) 

Since  an  integration  in  the  y-domain  comprises  three  well-defined 
regions,    one  should  separate  the  y-integral  into  three  complementary 
integrals,    one  per  each  region. 

Therefore,    analyzing  each  term  in  the  average  power  expression 
of  equation  {G12),    one  obtains, 


(146) 


(148) 
(149) 


+  '^|B"^V«)llc^'^'(oc)UmK)'2.Y  ^<ysh  Yzy  (150) 

(151) 

/  I  (153) 

4-2B^^'c*^)  C^%)  siob  Viy  <^^l  ^Y  ] 
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+  (B'"M  |B"'V«i1  h-C^'Ck)  Ic^Vrtl)  siih  Yju  a!5^,  ViW ] 


j^V'-fl-XA'-'MlA'^.^-k-^^'"/-"^ 


j^$fM)  =  -riD'"W|D<V)l€=M 


(154) 
(155) 
(156) 

(157) 

(158) 

(159) 
(160) 
(161) 

(162) 
(163) 
(164) 

(165) 
(166) 
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(167) 


'^  J168) 

Therefore,    the  expression  for  average  power  in  equation  (G12) 
becomes, 


and,    for  the  case  v/here  Y     is  imaginary,    the  following  relations  hold, 


'»  =  jY/ 


sinh  Yi(7  aj  si'o  Y^Y 

UTS  Ir^  V^  L/  =  C<JS  Y^    t/ 


1 


f 


>inh  'iY,cy  ritin  IV, 


c|=j..n  zr,y    ^ 
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Therefore,    equation  (170)  becomes, 


\coj^  (^'-y^^Y^  Id^^'c«)Ph-^o  W^hK)?)  e'^'^l 

-y  V  (|B.^^'(«)jC<%)-B'''c^)ic(^V^)|)  (170a) 

These  two  expressions  correspond  to  equations  (G13)  and  (G13a), 
respectively. 

The  y-dependence  of  the  average  power  disappears  through 
corresponding  integrations  on  regions  1,  2,  and  3,  as  indicated 
below: 

-  Region  1, 

-  Region  2, 
]  t,lf\W-Yi^d^  ^~Ls,inh-Z.Y^J^-'~di  (172) 


HYz 

j^cCiQY^^dLu.^  6;n^)2n,d  +  ic(  (173) 

j    s/ob  2VvwdY«^  (Ctfs^  2nci-|)  (174) 

/''sio"  Y^'cfA^.ia^  ^-,  sin  2Yz'ci  ( 175 ) 

fc^''Y^,^cjdUi.-la^-^^„s>in  2Y;'ci  (176) 
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,<!  (177) 

Region  3, 


jVMa.  = 


(178) 


-•o 


Substituting  equations   (171)  through  (178)  into  equations  (170)  and 
(170a),    one  obtains  equations  (G14)  and  (Gl4a). 
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/iPPENDIX    H    -    CHARACTERISTIC    IMPEDANCE    PROGRAM 
C 

c 

C  THIS  PROGRAM  SOLVES  THE  CHARACTERISTIC  IMPEDANCE  CF 

C  CCPLANAR  PARALLEL  STRIPS  ON  A  DIELECTRIC  SLBSTRATE. 

C 

C  PROGRAM    ASSUMES    FOLLOWING    PARAMETERS    ARE     HNCWK 

C  -   THICKNESS    OF    THE    SLAB, DD, IN    MILLIMETERS 

C  -    WIDTH    OF    THE     STRIPStWWtIN    MILLIMETERS 

C  -    DIELECTRIC'S    RELATIVE     PEP  MI TT I  VI TY t E R 

C  -    RATIO    OF    SEPARATION    BETWEEN    CONDUCTORS    TO    WIDTH    OF 

C  CONDUCTORS, SOW 

C  -    RATIO    OF    SUBSTRATE'S    THICKNESS    TO    FREE-SPACE 

C  WAVELENGTH, DRATIQ 

C  -    RATIO    OF    EFFECTIVE    DIELECTRIC'S    WAVELENGTH    TO    FREE 

C  SPACE    WAVELENGThtLRATIO 

C 

C  THE  METHOD  USED  IN  FINDING  THE  CHARACTERISTIC 

C  IMPEDANCE  IS  A  DIRECT  SUMMATiON  CF  INFINITESIMAL 

C  VALUES  BY  A  MODIFIED  SIMPSON'S  RULE. 

C 

C  THE    SUMMATION     IN     ALFA     IS    AN     APPROXIMATION    TC    AN 

C  INTEGRATION    IN    THE    ALFA-DOMAIN    FRCM     (-INFINITY)    TO 

C  (  +  INFINITY). 

C 

C  THE  APPROXIMATIONS  USED  IN  THE  PRESENT  PROGRAM  ARE 

C  -  THE  INTEGRATION  STEPtA 

C  -  THE  LIMITS  OF  THE  ALFA-DOMAIN  INT EGR AT  I  ON » AL FA 

C  AND  B. 

C  IT  IS  CLEAR  THAT  ANY  OF  THESE  VALUES  MAY  BE  CHANGED 

C  ACCORDING  TO  THE  SPECIFIC  CASE  OR  DESIRED  ACCURACY. 

C 

C  THE    QUANTITY    EO    IS    THE    FREE-SPACE    PERMITTIVITY    IN    MKS 

C  UNITS. 

C  THE    QUANTITY    MU    IS    th E    FREE-SPACE    PERMEABILITY    IN    MKS 

C  UNITS. 

C  THE    QUANTITY    C    IS    THE    FREE-SPACE    SPEED    CF    LIGHT     IN 

C  METERS    PER    SECOND. 

C  THE    QUANTITY    Ml     IS    THE    STRIP    CURRENT    IN    AMPERES. 

C 

C  THE  NUMBER  JJ  SETS  THE  NUMBER  OF  DRATIO  VALUES. 

C  ThE  NUMBER  OF  IL  SETS  THE  NUMBER  OF  SOW  VALUES. 

C 

C  PROGRAM    DEVELOPED    BY    LT(JG)     ARMANDO    LUNA    ECHEANDIA, 

C  PERUVIAN    NAVY. 

C  SUPERVISOR    -     PROF.     JEFFREY     B.     KNCRRtPH.Do 

C  U.S.    NAVAL    POSTGRADUATE    SCHOOL , MCNTE Re Y , C A L I F CRNI A 

C  SEPTEMBER, 1973. 

C 

c 

IMPLICIT    REAL=*^4    (M,N,K,L) 

DATA  WW/3.17/,DD/3.17/,ER/20./ 

DATA  A/0.5/ 

DATA    B/1700./ 

IL  =  7 

JJ  =  6 

PI=3. 14159 

FPI=4.vPI 

EC=l.E-9/( 36.*PI  ) 

MU=4.E-7' PI 

C=3.E8 

FACT0R=1.E9 

W=WW^1. E-3 

D=DD*l.£-3 

M1  =  W^^*2 

DC    11    1=1  ,IL 

READ    (5,101)     SOW 

101  FORMAT     (F5.3J 
S=SCW-W 

DC    13    J=1,JJ 

READ    (5,102)     DRAT  in,LRiTIO 

102  FORMAT     (F6.4,F8.6) 


FREQ=DRAT  lO^C/D 

FRE=FREC/FACTGR 

LAMBDA=C/FREQ 

LPRIME=LRAT1J^LAMBDA 

0NEGA=2.1  PI^  FREQ 

OE=OMEGA---E0-  ER 

CN  =  CMEGA--MU 

Kl=CPi£GA-SwKT(MU^EOJ 

K2=K1^SQRT{ER} 

K3  =  K1 

Kll=Kl*--2 

K22=K2-^2 

22    eETA=2.^PI/LPRIME 

B2  =  BETA-^-«2 

KCLSQ=K11-B2 

KC2SQ=K22-B2 

KC3SQ=KC1SQ 

KC4=(KC2SQ/KC3SQ) 

KC5=KC4-1 .0 

KC6=KC5**2 

M3=0.0 

ALFA=-1700. 

1    ALPHA=ALFA+^ /2. 0 

c 
c 

A2  =  ALPHA=^  t-2 

PHYSICAL    COMFIGURATI 

c 
c 

FCURIER    TRANS  FOR '^    GF 

c 
c 
c 


c 
c 
c 


IN    DEPENDEMCc     AS    DERIVED    BY    THE 
THE    CURREIVT    DENSITIES. 


Fl  =  ALPHA=MS+W)/2.0 

F2=ALPHA>W/2.0 

S1=SIN(F1) 

S2=SIN( F2) 

CJ=(4,/ALPHA)-.S1^S2 

AB=ALPHA- BETA 

G1=SQRT(A2+B2-K11) 

G3  =  G1 

ALGI=A2/G1-G1 

Gll  =  Gl^-i2 

G22=A2+e2-K22 

IF    (G22.LT.0.0)     GO    TO    7 

HYPERBOLIC    CASE,     GAiMMA    2     IS    REAL 

G2=SQRT (G22) 

G2D=G2-D 

G22D=2.-'G2D 

ALG2=A2/G2-G2 

ALAG2=A2+G22- 

SI=SINh(G2D) 

CC=C0SH(G2D) 

SI2=SINH(G22D) 

FCO  =  CCSF(G22D  )- 1 . 

INTERMEDIATE    ALGEBRAIC    STEPS 


NI  = 

N2  = 
N3  = 
N4  = 

14^C 
N'5  = 
Pl  = 
P2  = 
P3  = 
P^  = 
Ql  = 

14'^  C 
Q2  = 
Q3  = 

I*SI 
Q^  = 

l^SI 


KC4^^  ( 
KC4v  A 
AB/  (  0 
I./(C 

Oi 

iMI^-N4 

N4/N5 

N2/N5 

N'3/N5 

N1/N5 

KC2SQ 

0+KCI 

KC2SQ 

(AB-( 

-(  (OE 

((  AB^ 

+  (AB-"- 


G3/(G2-;  ER)--KC4  'SI+CO) 

&/(0E'^G2)^KC5-SI 

.M-^Gl  )^^KC^-^G3/(G2  ^ER)  "KC5--  SI 

MfGl)  *  {  (  Oy  >G2-(AB  .*^2}-  YsQ.b/  (CE- 

+  N2-N3 


G2)  J-SI+OM-'GO--KC 


-'  (  (  G3 /G2  )  *KC 4-='P  4- AB  /  {  OM  r  G  2  )^-  P 2'^>  KC  5 )  > S I W-KC2  SQ - P 

sg 

'  (  (AB/(0M-^G2)-^-KC5^  Pl-e(G3/G2)-^KC4'  F3  ) -S  I+P3^^  CC  ) 
G3/  G2  ) '  KC  4  < P4-UE  •  G 2>P  2-  {  A B-^ '<  2  )  /  (  OM  -  G2  )  -  KC  5  t^  P2 ) 
/ER)-  G3S^KC4-.  F2-AB-'P4-AR  -P4^=KC5  )  .^CO  +  AB 
^■2)  /(GM^-G2)  ^KC5  r^i  +  AB  M  G3/G2)  ^^  P3  ■■  KC4  +  GE-G2--P1 ) 
P3+(0E/Ek)^G3^KC4'P  L  +  Ab ;KC5-p3) vCC+ ( GE /ER )- G 1 
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c 
c 
c 


c 
c 
c 


Q5  = 
AE  = 
AH  = 
Cc  = 
Ci-  = 
BE  = 
Bh  = 
Dt  = 
DH  = 
AE2 
AI-2 
BE2 
BI-2 
CE2 
CH2 
Dt2 
CI-2 
Ti  = 
T2= 
T3  = 
TA  = 
T5= 
T6  = 
T7  = 
T8= 
T9  = 
TIC 
Til 
P=C 
1^T8 
GC 


g2*Q 

Ql/Q 

Q2/W 

Pl'i-A 

P3-*^A 

!./( 

l./( 

KC4V 

KC4= 

=  AE--^ 

=  AH-^ 

=BE 

=BHv- 

=  CE-^^ 

=CH- 

=DE^- 

=  DH*^ 

(AE2 

(AH2 

(BE2 

(6H2 

(CE2 

(CH2 

BE-^C 

B£>.-C 

Tl-M 

=  T3+ 

=T5  + 

MEGA 

^ALG 

TO    9 


3-Ql'' 
5*CJ 

E  +  P2- 
E-P4'- 
CE-*^G2 

0M^--G2 

CE 

CH 

>i  2 

■^z 

-2 
^2 
^2 

^Z 
+  DE2) 

+DH2) 

+  CE2) 

+  CH2i 

-BE2) 

-BH2) 

H-BH->' 

E*EO^ 

2 

T4 

T6 

■^BETA 

2)+2. 


Q4 


AH 

AH 

)»■  (  (CE/ER)  :<G3  /KC4v  CE  +  AB^KC5-^CI-) 

)^MABtKC5*CE  +  0M--G?^^KC4-t-CH) 


>MU 
-EG*  ER 

■^  EC*ER 

CE 
ER+MU*BH.-CH 


:  (ALG1>T9+0.5>^ALG2^  S  I2-*T10  +  D=»-T1 1-:  ALAG2  +  FC0 
iD> ALPHA 'G2>T7- (32-K2  2) 


TRIGONOMETRIC    CASE,    GAKiMA    2     IS    INAGiNARY 

G22=ABS(A2  +  B2-K2?.) 

G2=SQkT(G22) 

AG2=A2/G2+G2 

ALGG2=A2-G22 

G2C  =  G2-.  D 

G22D=2.-tG2D 

TSI=SIN(G2D} 

TCCJ  =  COS  (G20) 

TSI2=SIN(G22D) 

TC02=1.-C0S(G22D) 

INTERMEDIATE    ALGAERAIC    STEPS 


Nl  =  KC^i 
N2=KC4 
N3=AB/ 
N4=l-/ 

1KC4-:  TC 
N5  =  N1- 
P1=K4/ 
P2=N2/ 
P  3  =  1^3/ 
P4  =  M/ 
Q1  =  KC2 

lP4v'TC0 
C2=KC2 

10) 
Q3=(AB 

I'-TSI-C 
Q^=( (A 

I'i'TSI+l 

11 
Q5=-Q.Z^- 
AE=CI/ 
Ah=g2/ 
CE  =  P1^^ 
Ch  =  P3=« 
BE=1./ 


(G5/{G2  -ER) <KC4-TSI  +  TC0) 
AB/(Gt-  G2  )-;KC5'  TSI 

^GIJ  .KC4*  G3/(G2> -R)    KG  5-  TSI 

^"  G 1 )  -V-  (  (  -CM  r-G2-  (  A  D-'  -^2  )  -^KC 6/  (  GE  -G  2  )  )  -  T  SI  i-G  M-  G3^ 


{  GM 

(Ch 

0) 

N4+N2-''  N3 

N5 

fM5 

N5 

N5 

SQ- 

+  KG 

SQ- 


(  (G3/G2J  ^  KCAt  P4-/^B/(uM>G2)^:  P2^  KC5  )^  TS  H-KC2SG-' 

ISQ 

(  (<i&/(G.^i?-G2)*KC5*Pl  +  (G3/G2)    KC4-*F3  }  -TSI+P3-TG 


•-   ( G3  /  G 2  )  ••-  K C  4    P  4  +0  E  ■•  G 2vP  2-  (  A  B->  ^  2  )/(  0  M^  G2  )  -^^  KG  5^  P  2 ) 

(  CE/ER)-' G3^KG4.  F2-A6^  P4-aE-   P4-  KG5  )  fTCC  +  AB 

B-t--2)  /(0MvG2)'^KC  ?•;  PH-A5.  (  G?/G2]  -  P3^' KC4-GE^  G2<- PI) 

ae^  p3+  (  oe/  e  p  }  ^  g3  -  kg  4  .-  p 1  -j-a  &:  kcs^- p3  )  v  tg  c-*-  (  oe  /  er  i  -g 
Q3-qi-;q-4 

Q5-  GJ 

Q  5-  G  J 

AE  +  P2-^AH 

AE-P4-  AH 

(  0E^^G2i'M  {  CE/ER)^-G3*KC4^  GE  +  AE-.  KCS^^CI-) 
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BH  =  1./  (0M-1«G2  )*(AB»-KC5»  CE+OM^- G3  vKC4 

DE=KC4tCE 

DH  =  KC4-^CH 

/iE2=AE-^2 

Ah2=AH->T-2 

BE2=BE^*2 

EI-2  =  BH^^>-2 

Ct2=CE'^^*2 

DE2=DE:^--i2 

DE2=DH*-*2 

T1=(AE2+DE2)  --EO 

T2=(AH2+DH2)*MU 

T3  =  (B£2+CE2)-^E0<ER 

T4=(BH2+CH2)*MU 

T5  =  (CE2-BE2r'-^E0-^ER 

T6=(CH2'-BH2)^MU 

T7=BE^-CH-BK  CE 

T8  =  BE-^CE-E0-^ER+MU-BH^  CH 

T9=T1+T2 


CH) 


T10=T3+T4 
T11=T5+T6 
P  =  GMEGAiBETA-r< 
1 0^-D' ALP  HA-- G2-t- 

S    M2=A6S( P/FPI ) 
M3=M3+M2-A 
ALFA=ALFA+A 
IF     (ALFA.GT.B) 
GO    TO     1 

2    M4=y3/Ml 

WRITE     (6,19)     ER 


(ALGl'  T9  +  G.5^AG2.T10-*TSI2-vD^ 
T7  ••(B2  +  K22) 


GO    TO    2 


ALCG2^>Tll)-2 


19     FORMAT     ( 


3X, 


SCW,DRATI0,FRE,f^4 
FOR    ER    =     ',f5 


I  c 


tF5.2,2XT 'D 


1/L    =    •  ,F6.4,2Xt'  FREQ    =     ' , F5 .2 , 1 X t *  GHZ •  , 2X , • Z     IS     f,E14. 
17, IX, 'OnKS',//) 


13    CONTINUE 
11    CONTINUE 
23    STOP 
END 
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